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...Abstract. The total energy E of a star as a function of its angular momentum J and 
mass M in the Newtonian theory: E = E(J, M) [in general relativity, the gravitational mass 
M of a star as a function of its angular momentum J and rest mass m, M = M(J, m)], is 
used to determine the remaining parameters (angular velocity, equatorial radius, chemical 
potential, etc.) in the case of rigid rotation. Expressions are derived for the energy release 
during accretion onto a cool (with constant entropy), rapidly rotating neutron star (NS) in 
the Newtonian theory and in general relativity. A separate analysis is performed for the 
cases where the NS equatorial radius is larger and smaller than the radius of the marginally 
stable orbit in the disk plane. An approximate formula is proposed for the NS equatorial 
radius for an arbitrary equation of state, which matches the exact one at J = 0. 
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Introduction 



The change in mass and angular momentum of a cool neutron star (NS) during accretion leads 
to a transformation of its equilibrium state. Below, we consider NSs with weak magnetic 
fields that do not affect the accretion dynamics; i.e., our results refer to the accretion pattern 
in low-mass X-ray binaries and, in particular, in X-ray bursters. Of crucial importance is the 
question of what part of the energy released during accretion dissipates near and outside the 
stellar surface (in the accretion disk, in the boundary layer, in the spread layer, in the settling 
zone as matter is compressed under the weight of the newly supplied material, when uniform 
rotation is established by viscous forces throughout the entire extended stellar atmosphere, 
in the zone with a surface density of the order of 10 9 g/cm 2 subject to nuclear burning during 
X-ray flares) and what part dissipates inside the star in its interior. Observationally, these 
two zones of energy release differ radically. The energy released near the surface leaves the 
surface layers (is emitted) in a very short time: from fractions of a millisecond to several 
tens of seconds, although the complete energy dissipation can last for several hours. By 
contrast, the energy release in the stellar interior produces radiation with characteristic 
times exceeding hundreds of thousands and, possibly, millions of years; this is the time it 
takes for the interiors of neutron stars with weak magnetic fields to cool down (Levenfish et 
al. 1999). Such radiation can be detected only in highly variable transients when accretion 
on them virtually ceases. 

Here, we disregard energy release in the stellar interior by assuming that the star is cool 
and that its entropy does not change during accretion and transformation of the NS internal 
structure. This strong assumption allowed us to obtain a number of general results, which 
we used previously (Sibgatullin and Sunyaev 2000). Note that this assumption does not hold 
in the so-called thermal neutron stars, where the spinup via accretion leads to a difference 
between the angular velocities of the crust and the central liquid superfluid core inside the 
star. As a result, energy is released inside the star through viscous friction [see Alpar (1999) 
for a discussion]. 

A thermodynamic relation between the change in NS total energy (gravitational mass in 
general relativity) and the change in its angular momentum and mass (rest mass in general 
relativity) is considered for cool, rigidly rotating NSs with a given equation of state (EOS). 
This relation and the energy conservation law are used to derive a formula for the energy 
release on the stellar surface (in the Newtonian approximation and in general relativity) 
during disk accretion onto a NS rotating with an arbitrary angular velocity. The case where 
the NS equatorial radius is smaller than the radius of the marginally stable orbit is also 
considered in terms of general relativity. 

For the equilibrium figure of a rotating, incompressible fluid in its own gravitational field 
(Maclaurin spheroid), we derive explicit formulas for the dependence of disk and surface 
energy release on spheroid eccentricity. The fraction of disk energy release in the total 
energy release is shown to be expressed by a simple linear dependence on the ratio of the 
Maclaurin-spheroid rotation frequency to the Keplerian equatorial particle velocity. 

A universal approximate formula for the NS equatorial radius is derived for an arbitrary 
EOS. 

Previously (Sibgatullin and Sunyaev 2000), we considered astrophysical implications of 
our results and methods for deriving simple approximation formulas. We also analyzed the 
universal geometric properties of space-time outside rotating bodies. 



1 Newtonian treatment 

In the Newtonian approximation, the total energy of a NS consists of the gravitational, 
kinetic, and internal energies: 

E = J + \n\x 2 + y 2 ) + u{p{r)))p{f) dV (1) 

where the gravitational potential $ is 

$( f ) = G [ , p(r/) , dV (2) 
Jv\r-r/\ 

We take a model of an ideal gas with constant (zero) entropy for the internal energy of 
superdense matter. The Gibbs identity then yields 

du = —pd(l/ p) + Tds = —pd(l/ p). (3) 

We have the following obvious expressions for the NS mass and angular momentum: 

M= J p{f)dV, J = [ p(r)(x 2 + y 2 )ttdV (4) 
Jv Jv 

The following integral holds in a steady equilibrium state: 

u{p) + - - Vi 2 [x 2 + y 2 )/2 -$ = // = const, (5) 

since p has the same value at any point of the star [J because the equilibrium is isentropic, 
and because the rotation is rigid [the theorem of Crocco (1937); see also Oswatitsch (1976), 
Chernyi (1988)]. So, the dynamical equilibrium conditions can also be written as Vp = 
[see the case of an incompressible fluid in Lamb (1947)]. 

Let us consider two close equilibrium states with global parameters E, M, J, Q and E + 
<5E, M + 5M,J + 5J,tt + 6n. 

Theorem 1. The variations of NS mass 5M , angular momentum SJ, and total energy 
<5E for two close equilibrium states are related by the thermodynamic relation 

5E = n5 J + p5M (6) 

Indeed, denote the local displacement of point r on the stellar surface along the normal 
when passing from one equilibrium state to the other by W(f). For an arbitrary integral 
A = J v a(f) dV, we then have 

5A= [ 5a(f)dV+ [ a(f)W(f)dS (7) 

JV JdV 

In view of Eq. (2) for the gravitational potential, the variation of gravitational energy 5E gr 
is 

SE gr = -i / 5p$dV - ~ / p5$dV - \ f P W$dS. 



2 Jv 2 J v 2 jQ V 



11 The constant /i should not be confused with the Bernoulli integral which is constant along a 
streamline: [i = iq — [r, S7] 2 ! 



Reversing the order of integration in the second term on the right-hand side, we reduce this 
expression to 

SE gr = - [ 5p$dV- [ pW<$> dS (8) 

JV JdV 

We calculate the variation of internal energy SE in from formula (7) by using the Gibbs 
identity (3): 

SE in = 5 [ u(p)pdV = [ (u(p) + ^-)5pdV + / pW(u(p) + dS (9) 
Jv Jv P Jav P 

In order to simplify the subsequent calculations, we added the surface term J av WpdS to 
the right-hand part of (9); this term is zero, because the pressure vanishes at the stellar 
boundary. Clearly, the variation of kinetic energy can be expressed as 

5E = l I (5ptt + 25ttp)(x 2 + y 2 ){ldV +1 [ tl 2 p(r)(x 2 + y 2 )W dS. 
2 Jv 2 Jqy 

As a result, the variation of total energy 5E = 5~E gr + 5E in + <5E C is 

5E = J y {5p (-$ + ±fi V + y 2 ) + u(p) + p/p) + 5m(x 2 + y 2 ) } dV + 
f dv pW (-$ + u{p) + p/p + ln 2 (x 2 + y 2 )) dS. 

SM, 5 J are given, respectively, by 

5M= [ 5p(f)dV+ [ p(f)W(r)dS; 
Jv Jav 



5J= [ (n5p(f) + 5Qp(r))(x 2 + y 2 )dV+ [ Qp(f)(x 2 + y 2 )W dS. (10) 

Jv JdV 

Regrouping the terms, we represent the expression for <5E as 

5E = J v {5p (-$ - \n 2 (x 2 + y 2 ) + u(p) + h) } dV + 
j v {5pQ + 5Qp} {x 2 + y 2 )QdV + j dv Qp{r){x 2 + y 2 )W dS + 

! dv pW{-^ + u{p) + ^-\n 2 {x 2 + y 2 )) dS. (11) 

The coefficients in front of Sp in the first and last integrals of Eq. (11) match the constant 
p given by (5). Let us now make use of Crocco's theorem and factor this constant outside 
the integral signs. The angular velocity f2 can also be factored outside the integral signs in 
the second and third integrals, because the rotation is rigid. Having done these operations 
and using formula (7) for 5M and 5 J, we obtain 

5E = Q5 J + p5M, (12) 

hence, 

^ <9E, <9E, . . 

V=qj\m, P= m \j- (12«) 



4 



We see that the constant /x has the physical meaning of chemical potential here (Landau 
and Lifshitz 1976). 

Theorem 1 is conceptually associated with the following remarkable variational principle 
of general relativity by Hartle and Sharp (1967) (see also Bardeen 1970): the true mass and 
angular-momentum distributions differ from their virtual (possible) distributions with a fixed 
rest mass and a fixed total angular momentum in that they give a conditional extremum to 
the gravitational mass (energy in the Newtonian approximation). In this case, the angular 
velocity fl and the constant /i act as the Lagrangian factors. 

Corollary 1. When a star loses its angular momentum and energy by the radiation of 
electromagnetic or gravitational waves, the rates of change in its total energy and angular 
momentum are related by E = QJ. 

This equality follows from (12), because the change in NS baryonic mass is zero during 
wave emission. Corollary 1 was proved by Ostriker and Gunn (1969) when considering the 
fluxes of angular momentum and energy of electromagnetic or gravitational radiation in the 
wave zone. It has important applications for radio pulsars: the quasi-equilibrium evolution 
of the NS structure when it loses its angular momentum does not lead to any heating of the 
stellar matter. 

Denote the fluxes of angular momentum and energy on the NS by Ml and Me, respec- 
tively, where / and e have the meaning of the specific angular momentum and specific energy 
brought by the accreting matter. 

Theorem 2. When the angular momentum Ml and energy Me, are transferred to a 
neutron star by accreting particles in unit time, the following energy is released in the star 
with constant entropy in unit time: 

L 8 = M(e-ni-n), (13) 

where \i is the chemical potential of the cool star. 
Indeed, according to theorem 1, we have 

dE = QdJ + fidM. 

On the other hand, it follows from the energy conservation law that 

dE = dMe - L s dt, 

where L s is the rate of energy release (stellar luminosity) during accretion. 
From the law of conservation of angular momentum, we have 

dJ = Midi. (14) 

By equating the expressions for dE from theorem 1 and from the energy conservation law 
and using (14) for dJ, we reach the conclusion of theorem 2. 

Note. For the constant /i, we may choose its value at the stellar equator: \x = — \VL 2 R 2 — 
$ e , where <£> e is the gravitational potential at the stellar equator, and R is the equatorial 
radius. Here, we make use of the fact that the enthalpy vanishes on the stellar surface. This 
choice of \x allows us to determine the equatorial radius, the most important NS parameter 
[see formulas (32) and (33)]. 

Corollary 1. When a thin accretion disk is adjacent to the stellar equator and when the 
star transforms into an equilibrium state with a new mass and angular momentum in time 
dt, the energy ^dM(Qx — £l) 2 R 2 is released on the stellar surface, where dM is the amount 

Ft 



of baryonic mass accreted onto the NS in time dt; VLk an d f2 are the Keplerian velocity at 
the NS equator and its angular velocity, respectively. 

Indeed, in this case, e = v 2 /2 — $ e ; v = RQk,1 = R 2 &k, and expression (13) for L s 
turns into a full square. In particular, it follows from corollary 1 that the energy release for 
counterrotation at = —0.5^ is a factor of 9 larger than the energy release for corotation 
at Q = 0.5tt K \ 

The local justification here is as follows. Falling on the stellar equator, a particle of 
mass m increases the NS moment of inertia. An additional work is done, which is equal 
to the difference between the particle angular momenta in the Keplerian orbit, mR 2 fl K , 
and on the stellar surface, mR 2 il, multiplied by the NS angular velocity fi. This work is 
equal, with an opposite sign, to the additional (rotational) energy that must be added to 
the difference between the particle kinetic energies in the Keplerian orbit, mR 2 Q 2 K /2, and 
on the stellar surface,mi? 2 f2 2 /2. Therefore, when a particle approaches the NS from a thin 
accretion disk and decelerates in a narrow boundary layer at the equator from the Keplerian 
angular velocity VL K to the NS rotation velocity f2, it releases the energy l/2m(Vl K — il) 2 R 2 , 
giving up its angular momentum and part of its energy to the star. The formula for the 
energy release in this form was justified by Kluzniak (1987). The various local derivations 
of this formula were discussed by Kley (1991), Popham and Narayan (1995), and Sibgatullin 
and Sunyaev (1998). 

It should be noted that, having accreted at the equator, the matter cannot remain there 
for long. It must spread in some way over the surface while changing its angular momentum 
and energy. Therein lies the inconsistency of the local approach, which disregards the subse- 
quent redistribution of accreting matter over the star. The close match between the energy 
release given by (13) and the energy released when particles decelerate from the velocity 
RVL K to the velocity RVl on the stellar surface is fairly unexpected. 

Here, we make an attempt to solve the problem by taking into account the restructuring 
of a cool star during accretion. This global approach leads to formula (13). Thus, (13) holds 
not only for a narrow boundary layer near the equator (Popham and Sunyaev 2000). Our 
derivation of corollary 1 from theorem 1 suggests that the same energy is released by a particle 
in the spread layer (Inogamov and Sunyaev 1999), although the situation in the spread layer 
differs radically from the problem with a thin equatorial boundary layer. In the surface layer, 
the matter heats up, the thermal energy radiates away in bright latitudinal rings, and other 
complex physical processes accompanying the spread of matter over an equipotential of a 
rapidly rotating star take place. However, this all affects the instantaneous equilibrium of a 
neutron star only through global changes in its baryonic mass and angular momentum. It is 
pertinent to recall here that, according to Huygens's theorem (see, e.g., Tassoul 1978), the 
work done by gravitational and centrifugal forces on a particle to move it over the surface 
of a rigidly rotating liquid body in its own gravitational field is zero. 

Corollary 2. Since the kinetic energy Mv 2 /2 of the decelerated matter is released during 
spherical accretion onto a nonrotating cool star at rate M in its surface layer in unit time, 
there is no volume energy release in this case either. 

This assertion follows from theorem 2, because, in this case, / = Q = 0, $ = $ e = GM/R, 
and the total energy release given by (13) matches the surface one. 

Disk Energy Release. Passing from one Keplerian orbit to another with a smaller 
radius, an accreting particle loses its energy through viscous friction, which radiates away. 
The total energy radiated away by one particle from an infinitely distant point before it 
approaches a Keplerian circular orbit on the stellar surface is equal to the particle energy 
in this orbit, with an opposite sign. Therefore, the following formula holds for the energy 



release of the entire disk: 

L d = -M(R 2 Q 2 K /2-<Z> e ). (15) 
Setting the sum of centrifugal and gravitational forces equal to zero yields 

The formula for the disk energy release follows from (15) and (16): 

I^Ii^U (17) 
2R dr 1 H v ; 

If there is no rotation, then L d = L s = MGM/2R (Shakura and Sunyaev 1973). 

In order to estimate the effect of NS rotation on the disk energy release, let us consider 
an example in which the stellar matter is modeled by an ideal, incompressible fluid. In this 
case, the axisymmetric equilibrium figure is a Maclaurin spheroid. The Dirichlet formula 
(see, e.g., Lamb 1947) can be used for the gravitational potential of a homogeneous ellipsoid, 
and the following expression can be derived for L d from (17) (below, e denotes the spheroid 
eccentricity, e = a/1 — c 2 /a 2 , c, a to be the smaller and larger semiaxes): 

•3GM/e 2 -l VT^\ 

= M ^T (— arcsm e + )■ (18) 

The Keplerian equatorial rotation frequency can be calculated from (16): 



f K = y/G^B(e); B{e) = \/ arcsin e - (19) 

The formula for the rotation frequency is (Lamb 1947) 

/ = y/Gpg(e)/27r; g{e) = Vl 7 - (3 - 2e 2 ) arcsin e - 3^-. (20) 

e J e 

According to corollary 1 of theorem 2 and formulas (19) and (20), the energy release on the 
NS surface is 



^ 2d2 „ V 3GM / , /arcsin e VT^e 2 ,/3-2e 2 . ^VT^ 



L s = -M(n K - WR' = M— ( y — — - y — — arcsin e - 3- g2 

(21) 

For the case of counterrotation one should change the sign before the second therm in round 
squares of expression (21). 

Using (18 - 21), we can obtain the following unexpectedly simple relation to estimate the 
fraction of radiation from the disk in the NS total radiation: 

L d /(L S + L d ) = 0.5(1 + ///*). (22) 

This formula is valid both for the positive values of / (the case of corotation) and for negative 
ones (the case of counterrotation). 
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Note that the eccentricity e — > 1 as / — > for an incompressible fluid. Therefore, 
the spheroid asymptotically transforms into a plane disk. However, large angular velocities 
cannot be reached: Maclaurin spheroids become unstable to quadrupole perturbations at 
e > 0.8127 (at / > 0.4326 ^Gp/2n) (Lamb 1947; Chandrasekhar 1973). 

According to (18) and (21), the total luminosity L s +L d in the stability range of Maclaurin 
spheroids satisfies the inequalities 

L s + L d ( M V 2/3 / p V 1/3 

0.0485 < s . — — '- * < 0.11 

Mc 2 \1.4M J Vl0 14 g/cm 3 / 

in the case of corotation and the inequalities 

L s + L d ( M \ - 2/3 ( p \ - 1/3 

0.11 < — — '- ~ < 0.1953 

Mc 2 \1.4M Q J Vl0 14 g/cm 3 / 

in the case of counterrotation. The speed of light (constant) was introduced here for con- 
venience of comparing the results in the Newtonian approximation and in general relativity 
(Sibgatullin and Sunyaev 2000). 

For the critical e = 0.8127, the ratio of angular velocities Q/Qk is 0.602. It thus follows 
from formula (22) (in the stability range of Maclaurin spheroids) for spheroid and disk 
corotation that 0.199 < L S /(L S + L d ) < 0.5. For a nonrotating star as was shown by Shakura 
and Sunyaev 1973 the ratio L s j (L s + L d ) is equal 0.5. The inequality 0.5 < L s j (L s + L d ) < 
0.801 holds for disk and spheroid counterrotation. 



2 Treatment in General Relativity. 

The metric of stationary, axisymmetric spaces invariant to the change t,<f>=> —t, —<p can be 
written in cylindrical r, z, <fi coordinates as 

ds 2 = a 2 dt 2 - b 2 (d(f) - udt) 2 - e 2a [dr 2 + dz 2 ), (23) 

where a, b, a are the sought-for functions of r and z. The gravitational-mass functional in 
the steady-state, axisymmetric case is (Hartle and Sharp 1967) 

Mc 2 = c 4 / — ^—^jd^drdz + c- 1 [ T o °v^rf0, drdz, (24) 

Joe 16 Tt(-* Jv 

where V and dV are the region occupied by the star and its boundary, respectively. The 
integral of the scalar curvature is extended to the entire three-dimensional space: 

~ . . (a 2 + uoVtb 2 ) , , 0n 

r °^ + V-(.-n)^ - p; V=5 = »»«». 

The NS nuclear overcompressed matter is an ideal isentropic gas with a given internal energy 
density, e = pu(p) , with de = (p+e)dp/ p. According to formula (24), the mass functional for a 
steady-state, axisymmetric star is determined by the functions p, a, b, and a, of independent 
variables r and z and by the constant fl. Let us consider two close equilibrium states A and 
A + 5A. Subtracting the mass functional in state A from the mass functional in state A + 5A 
yields 

5M = 2tt / 5p— + 5a— + 5a— + 5b— }drdz + 5Q— + c~ 3 / T°Wyf^dl. 
J V op oa oa y°° J oil J dv 



8 



The functions of the form ^ are variational derivatives of the functional M with respect to 
the corresponding function a(r,z). 

Since the metric coefficients and their first derivatives are continuous at the boundary 
of the region dV ', the surface terms from variations of the gravitational mass reduce to 
the integral over a distant surface at pseudo-Euclidean infinity, where they vanish. The 
variational derivatives of M with respect to the metric coefficients are zero by virtue of the 
Einstein equations. 

Using the equalities u°uo + vftu$ = 1, and = flu , one could derive the following 
Hartle-Sharp formula at fixed metric coefficients: 

5u° = - U4> (u°) 2 6Q. (25) 

If we vary the gravitational-mass functional (24) at given metric coefficients, then we obtain 

5Mc 2 = 2nc- 1 5 [ Tg^drdz = 2nc- 1 5 [ (e - T$Q)^/=g~dr dz. (26) 
Jv Jv 

We now use the expressions for rest mass m and angular momentum J (Hartle and Sharp 
1967) 

m = 2tt / pu°y/^gdrdz, 
Jv 

J = -2c-\ / T°y/^drdz 
Jv 

and the relativistic integral of dynamical equilibrium conditions 

p + e 
(pcu°) 

In view of (24), the following formula holds: 



jj, = const. 



5e = P±^5p = P -^r(5( P u°) - p5u°) = P -^5( P u°) + Tpn. 
p pu u pu" v 

The quantity 5M can then be easily represented after reducing similar terms in (26) as 

5Mc 2 = n5J + pSm (27) 

hence 

„ dMc 2 , dMc 2 , , , 

n = ^ri»» "=-as-i' (28) 

Theorem 3. For any two close isentropic equilibrium states of a NS, the variations in 
gravitational mass, angular momentum, and rest mass are related by (28). 

Energy Release during Disk Accretion onto a NS in General Relativity. Let us 
now use theorem 3 to calculate the NS radiation and evolution from equilibrium states. The 
constant in the relativistic case is also convenient to estimate from its value at the stellar 
equator: 

H = c/u° = c^/a 2 -(u- Q) 2 b 2 , (29) 
the metric coefficients (23) are taken at the stellar equator. 



In general relativity, the NS equatorial radius can be larger and smaller than the radius 
of the marginally stable orbit, which is equal to three gravitational radii in the absence of 
rotation. In the former case, we denote the energy and angular momentum of a test particle 
of unit mass in a Keplerian orbit with the NS radius by e and /, respectively, take into account 
the law of conservation of angular momentum dJ = Idm, and equate the expressions for the 
change in gravitational mass from theorem 3 and from the law of conservation of energy: 

(fll + n) dm = dm e — L s dt. 

We use the equalities e = c (%)o, / = c (uk)^, A* = c/u°, and = Qu°, and express L s 
as 

dm 

L s = M(e -Ql-fi)= /x(u • - 1)M, — = M. (30) 

(JjL 

Here, we introduced the designation • for a scalar product of the 4-velocity vector of a 
particle rotating with the NS angular velocity and the 4-velocity vector of a particle rotating 
in a Keplerian orbit at the same point on the stellar equator. Expression (30) for L s is always 
positive and becomes zero if the NS angular velocity reaches the Keplerian one at the stellar 
equator. In our previous paper (Sibgatullin and Sunyaev 2000), we gave approximation 
formulas for the dependences of M, e, /, and fx, on Kerr parameter and rest mass for 
particles in equatorial Keplerian orbits for a NS with EOS A and FPS. We used (28) to 
calculate the NS angular velocity and chemical potential. For stars with EOS A and FPS and 
with a gravitational mass of 1.4 M Q , we derived approximation formulas for the dependences 
of total luminosity L s + L d and L S /(L S + L d ) ratio on NS angular velocity by using (30). 

The case where the NS radius is smaller than the radius of the marginally stable orbit 
is characteristic of general relativity alone. In this case, the gravitational field does an 
additional work on the spiraling-in particles in the gap; no stable Keplerian orbits exist 
within the marginally stable orbit. Thus, the rate of energy release is 

L s = M(e* - - n) = fi(u ■ u* - 1)M, (31) 

where u* denotes the 4-velocity vector of the particle that fell on the stellar surface from a 
circular marginally stable orbit with the energy e* and angular momentum /* corresponding 
to this orbit. In paper Sibgatullin and Sunyaev (2000), we gave universal (valid for an 
arbitrary equation of state) approximation formulas for the dependences of e*, /*, R*, and 
fx on Kerr parameter and on dimensionless quadrupole coefficient b for the particles falling 
on the surface from circular marginally stable orbits. We constructed plots of total luminosity 
against angular velocity for a NS with gravitational mass M = 1.4M© by using (27) and 
(30). 

Note the fundamental difference between our formulas (30), (31) for the luminosity and 
formulas (16), (18) from Thampan and Datta (1998). These authors calculated the energy 
release produced by particles during accretion as a difference between the particle energy in 
the stable circular equatorial orbit nearest the body and the energy of the particles lying on 
the equator and corotating with the body. 

Using Theorem 3 to Calculate the NS Angular Velocity and Equatorial Ra- 
dius. As the NS parameters, we choose its rest mass and Kerr parameter j (dimensionless 
angular momentum cJ/GM 2 ). The methods of constructing the function M(j,m) were 
developed previously (Sibgatullin and Sunyaev 2000). 



in 



According to (28), the NS angular velocity at given rest mass m and Kerr parameter j 
can be determined from the known function M(j, m) by using the formula 



QGM/c 3 = - M,j [™ , . 

Another important formula that relates the constant \i to a derivative of the gravitational 
mass with respect to the rest mass at constant angular momentum follows from (28). Taking 
(29) for fi, we obtain 



MM, m \j 
M + 2jM,j 



\/a 2 -b 2 {Q- u) 2 , R = b, (32) 



where R is the geometric equatorial radius of the star (the equator length divided by 2@[pi]); 
the metric coefficients a and b are taken at the NS equator. 

In the static case, the NS radius can be determined from (32) by using the Schwarzschild 
metric (w = SI = 0, a 2 = c 2 (1 - ^): 

^ = - Ar- (33) 



GM 1 - M 



'rn 



Remarkably, formula (33) [if M(j, m) is substituted in it and differentiated at constant j] 
closely agrees with the numerical data from Cook et al. (1994) and Stergioulus (1998) for 
the equatorial radius of a rotating NS, to within a few percent:.?? ~ 2GM/(c 2 — (?M 2 m 
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